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£S) ■ The mean free path of neutrino in charged and neutral current reactions is cal- 

culated for inhomogeneous nuclear matter which is expected to appear in the crust 

o 

of neutron stars. The relevant cross section depends on Fermi and Gamow- Teller 
strength distributions, which are derived from the large-scale shell model calcula- 
tions within the self-consistent Skyrme-Hartree-Fock approach and in a relativistic 
mean-field model. The inhomogeneous nuclear matter is described in terms of cubic 
Wigner-Seitz cells, which allows for a microscopic description of the structures in 
the so-called pasta phase of nuclear configurations and provides a smooth transition 
to the limit of homogeneous matter. The influence of pasta phase, its microscopical 
structure and geometrical shapes on neutrino propagation is considered. 

d 

^ ! I. INTRODUCTION 

o 

The transport properties of neutrino play an essential role in the physics of supernovae 
core collapse and in the evolution of the newly born neutron stars. The most important 
ingredient of neutrino propagation calculations is the neutrino opacity in a wide range of 
densities. Both the charged current (CC) absorption and neutral current (NC) scattering 
reactions are important sources of the neutrino opacity. In earlier works on neutrino inter- 
actions with the homogeneous nuclear matter the noninteracting baryons were considered 
hi. Later the strong interaction was taking into account both in the non-relativistic and 
relativistic calculations (see, e.g., and ref. therein). It was shown that the 

neutrino opacities of interacting matter may significantly altered from those for the nonin- 
4j . However, the use of homogeneous matter is a good approximation, while 



teracting case 

the existence of the stable quasi-nuclei in the crust of neutron stars is energetically favorable 
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and must be taken into account. 

At low densities, the nuclei in matter are expected to form the Coulomb lattice embedded 
in the neutron sea, that minimizes the Coulomb repulsion between the protons. With in- 
crease of density the nuclear pasta structures occur and the stable nuclear shape may change 
form from spherical droplet to rod, slab, tube and bubble shapes js|. Roughly speaking, 
the favorable nuclear shape is determined by a balance between the surface and Coulomb 
energies. In the following under "pasta phase" we will assume quasi-nuclear structures with 
spherical as well as non-spherical shape, which are embedded in a neutron sea. 

Various attempts have been made to describe the ground-state structure of pasta phase 
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3,Q, 



based on Thomas- Fermi approximation [9(, Quantum Molecular Dynamics [10 . 
Hartree-Fock and Relativistic Mean-Field calculations (RMF) within the Wigner-Seitz (WS) 

nnn 

cell approximation [16|, |17|, [18( . Later the dynamical properties of pasta such as the response 
function and neutrino mean free path (NMFP) were investigated [lO,[l^,Q, 3|- ^ was f° un d 
that the coherent scattering of neutrinos on inhomogeneous matter significantly reduces the 
mean free path. The collective modes in neutron rich skin of pasta were also calculated in a 
random phase approximation by using spherical Skyrme-Hartree-Fock (SHF) method 19]. 
It was shown the appearance of the super-giant resonance mode at low excitation energy, 
which may affect the specific heat of the crust of neutron stars. 

The calculation of the NMFP in pasta phase presented in this work is based on Hartree- 
Fock calculations in cubic WS cell [if], [l7, [3], which allows for the description of non- 
spherical quasi-nuclear structures such as rods or slabs and contains the limit of homogeneous 
matter in a natural way. The self-consistent calculations are performed for /^-stable matter 
in a density range for which the quasi-nuclear structures discussed above are expected to 
appear. For the nuclear Hamiltonian we consider Skyrme forces (SLy4) but also perform 
calculations within the relativistic mean-field (Hartree) approximation. The stability of the 
pasta phase with increase of the temperature is also discussed. 

The NMFP is extracted from the relevant cross sections of neutrino on different pasta 
structures. We pay special attention to the dependence of our results on the internal struc- 
ture of the pasta phase and its geometrical shapes. The mean free paths obtained from these 
inhomogeneous structures are compared with those calculated for homogeneous nuclear mat- 
ter at the same global density, thus one can estimate the influence of the inhomogeneous 
phase on the propagation of neutrinos. 
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After this introduction the details of Skyrme-Hartree-Fock approach with pairing will be 
outlined and a method to calculate the NMFP will be reviewed in section 2. In section 3 
we discuss the density dependent relativistic mean-field (DDRMF) model. The numerical 
results are discussed in section 4 and the final section 5 contains the main conclusions. 



II. SKYRME-HARTREE-FOCK CALCULATIONS 



A. Energy functional 
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he Skyrme-Hartree-Fock approach has frequently been described in the literature 20, 
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23] . Therefore we will outline here only a few basic equations, which will define the 



nomenclature. The Skyrme model is defined in terms of an energy density 7~t{r), which can 



be split into various contributions 21 



24] 



Ti = T~Lk + H0 + H3 + 7i c s + Tifin + T~L so + T^coui, 



(1) 



where Hk is the kinetic energy term, Ho a zero range term, H3 a density dependent term, 
TC e s an effective mass term, Tt^ n a finite range term and 7i so a spin-orbit term. These terms 
are given by 

Hk = 2^' 

H = \t [(2 + x )p 2 -(2x + l)(pl + pl)], 
K 3 = ±t 3 p a [(2 + x 3 )p*-(2x 3 + l)(p 2 p + pl)}, 
Heg = l[t 1 (2 + x 1 ) + t 2 (2 + x 2 )]rp 

+ | [* 2 (2x 2 + 1) - t 1 {2x 1 + 1)] [r pPp + 
ftfin = -^[3ti(2 + xi) -t 2 (2 + x 2 )]pAp 

+^ [3t 1 (2x 1 + 1) + t 2 {2x 2 + 1)] [p p Ap p + p n Ap n ] , 
n so = -\W [ P VJ + Pp VJ p + Pn VJ n }. 



(2) 



The coefficients tj, Xi, W , and a are the parameters of a generalized Skyrme force 
imaginary time step was used to solve the Hartree-Fock equations 
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25|. The 



171 . |26j . The calcula- 



tions were performed in a cubic Wigner-Seitz cell with a size of typically 20 fm for charge 
neutral matter containing protons, neutrons and electrons in /3-equilibrium. Pairing corre- 
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FIG. 2: Profiles for the proton density distribution forming a slab-structure at a density of 

0.0775 fm- 3 . 

lations are included in terms of the BCS approximation by assuming a density-dependent 
zero-range pairing force, which has been used in earlier calculations [3, 27]. 

Figure [T] displays a typical example for the density profile of the proton distribution. At 
the global nucleon density of 0.0625 fm -3 we obtain a quasi-nuclear structure which can be 
characterized as rods along the z-axis. At a density of 0.0775 fm -3 the variational Hartree- 
Fock leads to a quasi-nuclear structure, which can be characterized as a set of parallel slabs. 
In figure 2, which shows the proton density distribution at this density, we have chosen the 
orientation of the coordinates such that these slabs are orthogonal to the z-axis. 

The single-particle energies and wave functions for protons and neutrons resulting from 
such Hartree-Fock calculations were used to evaluate the NMFP by using the method out- 
lined in the following subsection. 
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B. Neutrino mean free path in Skyrme-Hartree-Fock model 



The matrix element for the neutrino-nucleon reactions v + n 
is given by 



where 



_ GfC 



f = -m„y*(l - 7 5 )u v ( e ) 



v + n(v + n— > p + e ) 
(3) 

(4) 
(5) 



are hadronic and leptonic currents, respectively. The parameters A, V must be replaced 
by the respective values of coupling constants and C stands for the Cabbibo factor in the 



charged- current reaction 



. The total cross section can be written as 



M 



G 2 F C 2 

TV 



V 2 {1 + cos i?) I Mi 1 2 + A 2 {1 - - cos 



\M,\ 2 



(6) 

(7) 
(8) 



where within this non-relativistic approach we neglect the lower components in Dirac spinors 
it ~ (q). In the charged current reaction the influence of the Coulomb field on the outgoing 
electron can be taken into account by multiplying the cross section by the Fermi function 
F(Zf, E e ) 29|. For this reaction Mi and M 2 stand for the Fermi and Gamow- Teller matrix 
elements, respectively. The integration is performed over the spatial angle i9 between the 
momenta of incoming and the outgoing leptons. The single-particle wave functions <p{r) and 
single-particle energies £/ (ej) are obtained from the solution of the HF equations. Note 
that these single-particle energies enter ([6]) as the energy for the outgoing lepton is defined 
as 

E v ( e ) = E'J 1 + — 

where is the energy of the incoming neutrino. 

The formalism described so far is appropriate for the neutrino-nucleus interaction. With 
some extensions it may also be used to evaluate the interaction of neutrinos with the quasi- 
nuclear structures in the crust of neutron stars. Unlike spherical nuclei and the case of the 



6 



droplet phase, the cross section of neutrino on rods and slabs, generally speaking, depends 
on the spatial orientation of momentum transfer (fin (JSJ), since the density distributions of 
rod and slab phases are non spherical, as it is shown on figures fT|2l The precise averaging 
over all possible mutual orientations of vectors q and f requires additional numerical efforts. 
Thus, in order to reduce this effort we considered three particular cases, with the vector q 
along the direction of the x, y and z-axis. Doing so, we determine the averaged cross section 
as 

a = -(a x + a y + a z ), (9) 

where cr x (y,z) represents the cross section calculated for the momentum transfer along x(y, z)- 
axis. 

In contrast to a finite nucleus, the WS cell of the inhomogeneous nuclear matter contains a 
large number of unbound neutrons, which give nonzero contribution to the total cross section. 
Thus, the cross section consists of two parts: the cross section due to the interaction with the 
nucleons bound in the quasi-nuclear structure and the cross section due to the interaction 
with unbound neutrons. Therefore, one can consider (jf)]) as a cross section of neutrinos with 
all nucleons in a given volume V ce u of a WS cell. The reverse NMFP can then be written as 

la . . 

(10) 



A V ce ii 

Another important distinction of pasta structures in the crust of neutron stars from the 
finite, isolated nuclei consists in the existence of the electron sea in the volume of the WS 
cell. Therefore the nonzero chemical potential of electrons must be taken into account in 
the evaluation of charged current reactions by a blocking of final states for electrons with 
energies below the respective Fermi energy /i e . 



III. RELATIVISTIC MEAN-FIELD CALCULATIONS 

In order to test the sensitivity of the results on the underlying nuclear model and the 
choice of the NN interaction we also investigated the dynamical properties of inhomogeneous 
nuclear matter evaluated within a relativistic mean-field (Hartree) approximation by using 
a model of density-dependent meson-nucleon coupling constants. The parameterization of 
these constants has been fitted to reproduce the properties of the nucleon self-energy eval- 
uated in Dirac-Brueckner-Hartree-Fock (DBHF) calculations of asymmetric nuclear matter 
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)ut has also been adjusted to provide a good description for bulk properties of finite nuclei 



3. 
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33J. The density-dependent relativistic mean- field (DDRMF) approach has also 



been used to describe the properties of inhomogeneous nuclear matter in the crust of neutron 
stars kfil. 



A. Density dependent relativistic mean-field approach 

The relativistic mean-field (RMF) approach is an effective field theory of interacting 
mesons and nucleons. The Lagrangian density consists of three parts: the free baryon La- 
grangian density Cb, the free meson Lagrangian density Cm and the interaction Lagrangian 
density £ int : 

£ = £ B + £ M + £ int , (11) 

which take the explicit form 

£B=y(il^-M)y, 

2 Yl (l F (K)H» F (K) - m l A (^ A (K] 
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(12) 



Ant = - 9a^^a^> ~ 98^T^ S ^ 

- gjt nfpAfa* - gfir^Afa* 
_ e tt 7M I(l + T 3 )^ 7) tt, 

with the field strength tensor F( K ^ U = d^A^ — d u A^ fl for the vector mesons. In the above 
Lagrangian density the nucleon field consisting of Dirac-spinors in isospin space is denoted 
by \l/ and the nucleon rest mass by M = 938.9 MeV. The scalar meson fields are $ CT and 
$5, the vector meson fields Ai u -\ and At p y The mesons have rest masses m K for each meson 
k and couple to the nucleons with the strength of the coupling constants g K , which depend 
on a density of the nucleon field A very convenient parameterization for this density 



dependence has been given in 



32] 



The numerical procedure to solve the Dirac equation in the cubic WS cell is the same 



as m 



3, 0, Q 



Pairing correlations are included in terms of the BCS approximation 



assuming a density dependent zero-range pairing force, which has already been discussed. 
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The resulting single-particle energies and spinors were used in the calculation of NMFP as 
described in the next Subsection. 



B. Neutrino mean free path in relativistic mean- field model 

First, let us consider the charged current reaction. Here we will exploit the most general 
form for the nucleonic current, which is allowed due to the Lorentz, parity and isospin 



invariances 28J 



tf c = i^F^q 2 )^ + Fl{q 2 )a^q v + F A {q 2 )l^ - iF^q 2 )^}^ (13) 

where F% and F% are isovector electromagnetic formfactors, Fa is the axial-vector formfactor, 
F p is the induced pseudoscalar formfactor. Following the common practice we ignore the 
contribution of the second-class currents. The leptonic current has the same structure as in 
(JSJ). Analogously to (j3J), (JTj) the averaged squared matrix element for the charged current 
reaction can be written in the form 

W\ 2 = ^-[|-Mi| 2 (l-^rCos^) + |^ 2 | 2 (l + |^cos^) (14) 

+ \M z \ 2 {p 2 + E 2 V - 2 Pl E v costf - |L cos $ -^E 2 V cos + M^) 

+ \M 4 \ 2 ^((p 2 + El) cos $ - 2E l E v aaa& + -p l E l + —E 2 U - 2 Pl E v )], 
Ei pi 

where 

Mi = Fi4> p ji/; n + F A $ p j 5 ji/; n , 

M 2 = ^lV'pToV'n + FA^p^Qipn - iFp4> p ~f 5 q tp n , 

M 3 = F^ v ti\) n: 



M A = Fplfjp^lpn, 



and 



a 
a 

Dirac spinors ip and the respective single-particle energies are obtained from the solution 



of the Dirac equation 



3- 
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The hadronic part of the neutral current involves additionally isoscalar electromagnetic 
form factors F* and _F|, so that 

J»° = ^n[F A (q 2 h 5 j, - iF p (q 2 ) l5 q, (15) 
+ (1-2 sin 2 Ow) (F?(q 2 ) % + F v 2 {q 2 )a, v q u ) 
-2sin 2 6 w (F°(q 2 )^ + F^(q 2 )a^q v )}^ n , 

where 6\y is the Weinberg angle. The respective matrix element for this reaction looks like 

W\ 2 = if[ |M| 2 (l-icos^) + |AT 2 | 2 (l + cos^) (16) 

+ \M 3 \ 2 (E' U 2 + El - 2E u E' u cosd - \E'?cosd - ±£ 2 costf + \E V E' 2 ) 
+ |Ar 4 | 2 (^-^) 2 (l + cos^)], 

where 

M = ((1-2 sin 2 9 w )F v l - 2 sin 2 6 w F?)$M n + F A $ nl5 j, 

■M> = -F A $nl5lo^n - ((1 - 2 sin 2 O w )Ff - 2 sin 2 Q w Fl)^ n i\) n + %F P ^ n ^ n , 

H-i = ((1-2 sin 2 9 W )F% - 2 sin 2 9 w F^ n EiJj ni 

M A = F p ip n -f 5 ip n . 

Substituting (THI) and ( |T6l) in (jBJ) one obtains the mean free path of neutrinos in relativistic 
mean-field model for charged and neutral current reactions, respectively. 

IV. RESULTS AND DISCUSSIONS 

Before we start the discussion of the transport properties of neutrinos in the crust of 
neutron stars let us review some details of self-consistent Hartree-Fock and relativistic mean- 
field calculations. For all the results presented in this manuscript a temperature T — 1 MeV 
was considered. This temperature is high enough to take into account some effects of finite 
temperatures and low enough to maintain the pairing correlations and stable quasi-nuclear 
structures. If the temperature rises, the pasta phase structures become smoother and at some 
critical temperature they disappear. For the Skyrme-Hartree-Fock approach this critical 
temperature is around 5 MeV and 10 MeV for slab and rod structures, respectively, while 
the droplet structure disappears at a temperature higher than 15 MeV. Thus, the spherically 
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density p [fm ] 

FIG. 3: Proton abundance in the case of uniform matter (dashed line) and pasta phase 
(solid line). The symbols refer to specific calculations, whereas the lines have been added 
to guide the eye. The results for Skyrme-Hartree-Fock calculations are shown in the upper 
panel and the relativistic mean-field results in the lower one. The dashed arrows indicate 
typical densities leading to pasta structures of droplet, rod and slab shape. 

symmetric droplet phase will play the main role in different simulations containing the 
temperature evolution. Employing the relativistic approach the pasta phase structures turn 
out to be less stable and even the droplet structures disappear at a temperature of T = 10 
MeV. 

For temperatures below 1 MeV and global densities below 0.08 fm -3 the variational 
calculations yield structures with inhomogeneous density distributions (see figures Q] and 
|2]). Comparing the spectra of single-particle energies obtained for the homogeneous and 
inhomogeneous solutions one observes that the single-particle energies for the localized states 
are more deeply bound than the corresponding single-particle states for the homogeneous 
approach. In the /^-equilibrium all proton states are localized and therefore tend to have more 
attractive single-particle energies in the inhomogeneous as compared to the homogeneous 
density calculation. The variational calculations allowing for pasta structure yield larger 
proton fractions than obtained for the /^-equilibrium of homogeneous matter at the same 
global density. 
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FIG. 4: The neutrino mean free path (NMFP) calculated for the charged current reaction 
in case of rod and slab configurations demonstrate the dependence of the result on spatial 
orientation of the momentum transfer q. For these calculations we have employed results 
of the Skyrme HF approach and ignore the blocking of final electron states. 

This can be seen from inspecting figure O The upper panel of this figure contains results 
of the proton abundances for baryonic matter in /3-equilibrium resulted from non-relativistic 
Skyrme-Hartree-Fock calculations. The proton abundance of homogeneous matter is a mono- 
tonically increasing function of total density and it reaches the value of 4% at the density 
0.1/m~ 3 . Allowing for inhomogeneous matter distribution one obtains a significant increase 
of the proton fraction at densities below 0.03/m -3 , while in the density region from 0.03 
to 0.08/m~ 3 its value is almost constant around 3.2%. The lower panel of figure [3] displays 
the corresponding results derived from the relativistic mean field approach. This relativistic 
approach seems to provide a smaller symmetry energy at these low densities, which leads to 
smaller proton abundances in the inhomogeneous as well as the homogeneous solution. 

In figure 0] we want to demonstrate the dependence of the neutrino cross section for the 
charge current reaction on the spatial orientation of the momentum transfer q. This is 
displayed in terms of the corresponding neutrino mean free path, which has been calculated 
according to (fTU]) from a x (solid line) and a z (dashed line), respectively. Note that due to 
our choice of the coordinate system the results for a y are identical to those for a x for the 
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FIG. 5: Skyrme-Hartree-Fock calculations of NMFP for pasta phases (solid curves) and 
the respective results for homogeneous matter at the same global density (dashed curves). 
The results for the charged current reaction are shown in the left column, while the neutral 

current NMFP in the right one. 

rod as well as the slab structures. 

For the density p of 0.0625 fm -3 , which leads to a rod structure, we obtain results for 
the NMFP ranging 20 km for neutrinos with an energy of 10 MeV down to 30 cm for 
neutrinos with an energy of 100 MeV. For low-energy neutrinos the NMFP for reactions 
with a momentum transfer parallel to the x-axis is larger by a factor of 2 as compared to 
a momentum transfer parallel to the z-axis, difference which disappears for neutrinos with 
larger energies. This factor of 2 is non-negligible but small on the scale of variations for the 
NMFP as a function of the neutrino energy. Therefore the simple averaging procedure of 
([9]) seems to be adequate. 

Similar results are obtained for the slab configuration as can be seen from the lower panel 
of figure HI Note that the results for the NMFP are considerably smaller at low neutrino 
energies (by a factor of 10) and even at neutrino energies as large as 100 MeV smaller by a 
factor 2, although the ratio of the inverse densities is only about 1.2. 

The NMFP calculated in CC and NC reactions for homogeneous and inhomogeneous 
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FIG. 6: Results of NMFP due to NC (dashed lines) and CC (solid lines) reactions. The 
description of the inhomogeneous baryonic matter distributions results from the density 

dependent relativistic mean-field calculations. As examples we present results for the 
droplet phase displayed in the left panel at a density of 0.034 fm~ 3 and for the rod phase 

(right panel) at 0.055 fm~ 3 . 

matter distributions are shown in figure [51 First, let us compare the NMFP of homogeneous 
matter for both types of reactions. The main influence on NMFP's results from the available 
phase space for each reaction. In fact, the proton fraction of homogeneous matter does 
not exceed 1% for the densities considered here. Thus we have to consider a much larger 
blocking effect for the neutrons in the final states NC reactions than for the protons in the 
CC reactions. Therefore the cross section of CC absorption is larger than in NC scattering, 
and consequently, the mean free path is shorter, as it is shown by the red dashed lines in 
figure O Due to the small proton abundances in homogeneous matter the Pauli blocking 
factor of final electron states affects the result for the CC reaction only at very small neutrino 
energies E v < 10 MeV. 

Figure [5] also presents results for the NMFP of inhomogeneous matter for both types of 
currents. First of all, we should emphasize the larger influence of electron blocking factor 
on CC current reaction in the droplet phase. This is due to the larger proton abundances in 
the /3-equilibrium of the inhomogeneous matter. At a neutrino energy around E u ~ lOMeV 
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the mean free path of CC processes is longer in comparison with NC scattering, because in 
this region the Pauli blocking of electrons in CC reaction dominates over the differences in 
phase spaces of the baryonic states. If the energy of incoming neutrino E u rises the Pauli 
blocking drops exponentially and the ratio of the cross sections for CC and NC reactions 
is determined by the available phase space for the baryonic states as discussed above for 
the homogeneous matter calculation. This means that the NMFP of absorption due to CC 
becomes shorter than the respective result in NC scattering. At higher densities, where rods 
and slabs appear, the influence of Pauli blocking of electrons is partially compensated by the 
effects of the baryonic phase space. Therefore the NMFP of CC reaction remains shorter in 
comparison with NC reaction for all neutrinos with 10 < E v < 100. 

The same features are also observed in the comparison of NMFP due to the different 
currents for the models of inhomogeneous baryonic matter, which are based on the relativistic 
mean field calculations displayed in figure [6j 

The cross section for neutrino scattering in homogeneous matter increases with the bary- 
onic density in a non-linear way (see discussion above). Therefore one may expect that 
the mean free path in the inhomogeneous matter is shorter than the corresponding one for 
homogeneous of the same global density, since the scattering on the quasi-nuclear structures 
shall enhance the respective cross section. Nevertheless, the NMFP obtained for the charged 
current reaction, which is shown in left column of figure [5] demonstrates the opposite be- 
havior, specially at low densities, where the droplet phase occurs. The NMFP obtained 
from absorption in inhomogeneous matter is longer than the respective result derived from 
homogeneous matter calculations. 

In order to explain this effect one should consider difference in proton fractions of homo- 
geneous and inhomogeneous matter discussed in the beginning of this Section. At a typical 
density 0.0165/m~ 3 , where the droplet phase occurs the proton abundance in inhomoge- 
neous matter is significantly larger than the respective value obtained in the homogeneous 
matter. This difference in the proton fractions has two effects: first, the homogeneous matter 
contains less protons in comparison with the inhomogeneous one. Consequently, the number 
of unoccupied final proton states is larger and more transitions, which contribute to the total 
cross section, are possible. Secondly, the chemical potential of electrons compensating the 
charge of the protons in matter is lower in case of homogeneous matter and the respective 
Pauli blocking factor for the produced electrons is lower than those obtained for the inhomo- 
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geneous matter. This effect again modifies the cross section considerably at low E v . With 
increase of the energy of incoming neutrinos the Pauli blocking of electrons rapidly drops 
and more transitions become possible, so that the differences between homogeneous and 
inhomogeneous matter distributions are getting less significant and the respective NMFP's 
become closer one to another. At higher densities of matter, where the rod and slab phases 
occur, the difference in proton abundances are less important, therefore the resulted mean 
free paths are very similar and the effect of inhomogeneous structure becomes negligible. 

At the end we should notice that at neutrino energies less than 10 MeV (thermalized 
neutrinos) the NMFP's of homogeneous and inhomogeneous matter distributions calculated 
in CC reaction significantly exceed the typical neutron star radius. Therefore one can 
conclude that the charged current reaction is kinematically suppressed [34]. 

The results of neutral current reaction are shown on the right panel of figure [51 It is obvi- 
ous that the appearance of pasta phase in this case has no important influence on neutrino 
propagation, since this type of reaction does not depend on Pauli blocking of neutrino in final 
state (no trapped neutrinos). The only small difference in NMFP's of homogeneous matter 
and droplet phase may be explained by different values of matrix elements in (jSJ), since the 
s.p. wave functions of bound neutrons in droplet significantly differ from wave functions of 
homogeneous matter. However even this small effect becomes negligible if the global density 
increases and the density profiles become smoother and transition to the homogeneous phase 
approaches. 

A comparison of NMFP's of charged and neutral currents in case of pasta phase based on a 
relativistic mean-field model in a WS cell is displayed in figureEl It is worth mentioning that 
within the relativistic model we could not find any formation of slab structures. Therefore 
only results for droplet and rod structures are shown. Also, the global density, at which 
the droplet phase occurs in the relativistic mean-field model is two times larger than the 
respective density in the nonrelativistic model. The difference between proton fractions of 
homogeneous matter and pasta phase is not so significant. In fact, the values of proton 
abundance around p = 0.02 fm -3 , displayed in the lower panel of figure El are about 40% 
smaller than the corresponding values obtained in the Skyrme model (the upper panel). 
Therefore we omit the comparison between NMFP's of homogeneous and inhomogeneous 
matter however we compare the mean free paths of pasta phase for both types of reactions. 
One can see that at E v < 20MeV the behavior of CC curves is determined by the Pauli 
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blocking, while at higher energies the result becomes sensitive to the structure of phase space 
available for the reactions. Both charged and neutral current mean free paths decrease if 
the global density of matter rises. 

Summarizing we conclude that the NMFP is determined by two different factors. The 
first of them - the Pauli blocking effect of final electrons in CC reaction play the most 
important role at low neutrino energy and drops exponentially if the energy increases. The 
second factor is the difference in baryonic phase spaces of different reactions. The phase 
space of CC absorption is larger than in NC scattering, because the Fermi energy of final 
(proton) states is considerably lower than the neutron Fermi energy. 



V. SUMMARY AND CONCLUSION 



The aim of this study was to examine the role of the inhomogeneous baryonic density 
structures in the crust of neutron stars on the propagation of neutrinos. Our calculations 
of neutrino mean free paths (NMFP) are based on microscopic descriptions of the so-called 
"pasta structures" derived from 3D Hartree-Fock calculations with the SLy4 parameteriza- 
tion of the Skyrme potential as well as density-dependent relativistic mean-field calculations, 
which reproduce the empirical properties of normal nuclei with good accuracy. We find that 
the evaluated NMFP due to charged current reactions significantly depend on the structure 
of the pasta phase. This is mainly due to fact that the proton abundances derived from the 
/^-equilibrium in inhomogeneous matter are larger than the corresponding values determined 
for symmetric matter. The effects of inhomogeneous baryonic density distributions is less 
pronounced for the neutral current contribution. 

Recent studies show that the weakly bound neutrons may play an important role in 
formation of collective modes in the crust of neutron stars 19(. In our calculations the 
role of such collective features of neutrons has not yet been considered and an accurate 
calculation of nuclear response functions should be done in the future. 
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